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ABSTRACT 

Dark matter in the Universe consisting of macroscopic objects such as primordial black holes 
may cause gravitational lensing of distant objects. The magnification associated with lensing 
will lead to additional scatter in the received flux from standard candles, and too small an 
observed scatter could rule out compact dark matter entirely. In this letter, we show how the 
scatter in fluxes of distant Type la supernovae could be used to distinguish between mod- 
els with and without lensing by macroscopic dark matter. The proposed SNAP project, with 
- 2400 supernovae in the range 0.1 <J z ^ 1.7, should be able to identify models at 99.9% 
confidence, if systematic errors are controlled. Note that this test is independent of any evolu- 
tion of the mean supernova luminosity with redshift. The variances of the current Supernova 
Cosmology Project sample do not rule out compact lenses as dark matter: formally they favour 
such a population, but the significance is low, and removal of a single faint supernova from 
the sample reverses the conclusion. 



1 INTRODUCTION 



In recent years one of the classic tests of the geometry of the Uni- 
verse has undergone a resurgence of interest. Type la supernovae 
are thought to act as standard candles, and therefore their Hubble 
diagram may be used to constrain cosmological models. Dedicated 
discovery and follow-up programmes (Perlmutter et al. 1997, Perl- 
mutter et al. 1998, Schmidt et al. 1998) have established the Hubble 
diagram to z ~ 1, leading to constraints on Ho, Qm and Qa, most 
notably the requirement of a positive cosmological constant (Riess 
et al. 1999). It has long been recognised that the supernova Hub- 
ble diagram would be affected by gravitational lensing, if signifi- 
cant quantities of dark matter resided in the form of macroscopic 
compact objects (herafter MACHOs) such as black holes (Linder, 
Schneider & Wagoner 1988, Rauch 1991, Holz 1998, Metcalf & 
Silk 1999, Perlmutter et al. 1999, Seljak & Holz 1999, Weller & 
Albrecht 2000, Huterer & Turner 2000). Such a population has 
been proposed to account for the long-term variability of quasars 
(Hawkins 1993). 

These studies have shown that lensing can have a significant 
effect on the estimation of cosmological parameters, and that the 
distribution of supernova fluxes could be used to determine whether 
the dark matter was in MACHO form or not. Metcalf & Silk (1999) 
showed that, provided the cosmology was known, microscopic dark 
matter could be distinguished from MACHOs with relatively small 
numbers of supernovae. The reason for this sensitivity is that in a 
MACHO-dominated universe of reasonable density, most bunches 
of light rays do not undergo large magnifications, and the most 
probable flux received is close to that expected in empty-beam 
models (Dyer & Roeder 1974). This shifts the most probable Hub- 
ble diagram systematically, changing the estimates of the cosmo- 
logical parameters (Holz 1998). Given that one does not know 
the underlying cosmology, one should estimate simultaneously the 
background cosmology as well as the contribution to the matter 
density by MACHOs. This is ambitious, and will surely be at- 



tempted once larger supernova searches such a s those proposed by 
the Supernova/ Acceleration Probe (SNAP: see http://snap.lbl.gov) 
and by the Visible and Infrared Telescope for Astronomy (VISTA; 
http://www.vista.ac.uk) are underway. In this paper, the study is 
more limited; we focus on how the extra scatter in supernova fluxes, 
from lensing, can be used to constrain the quantity of MACHO dark 
matter. The variance will have contributions from intrinsic varia- 
tions in supernova properties, instrumental error, and lensing. The 
first two of these should be virtually independent of the cosmolog- 
ical model; the presence or absence of significant additional scatter 
allows a test of whether the dark matter is in the form of MACHOs. 

The detectability of dark matter candidates in this test rests 
on how accurately the variance can be estimated from sets of su- 
pernovae. For gaussian distributions, this is readily done, but, al- 
though it may be a reasonable approximation to model the intrinsic 
and instrumental effects as gaussian, the lensing effect is highly 
skewed towards rare high-magnification events. We must therefore 
take care in modelling this effect accurately. We use numerical ray- 
tracing simulations (described in section 2) and analyse the results 
using a Bayesian method (section 2.1). Lensing induces an ex- 
tra scatter which rises with redshift, contributing as much as 0.5 
magnitudes at a redshift 1.5 for an Einstein-de Sitter model with 
all the matter in MACHOs. This is readily distinguishable from 
a no-lensing model, provided sufficient supernovae are available. 
We find, in general, that with ~ 2400 supernovae selected from 
0.1 < z < 1.7, as proposed by SNAP, the scatter alone can select 
between several cosmological models at 99.9% confidence level, if 
systematic errors are controlled. Note that this test has the advan- 
tage that it is insensitive to any evolution in the mean luminosity of 
supernova (although it is dependent on any evolution in the scatter 
of intrinsic properties). Similar conclusions to ours have recently 
been reached via a slightly different analysis by (Mortsell, Goobar 
& Bergstrom 2001). 
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Figure 1. The ray-tracing simulation geometi'y. In practice we use 5 planes, 
but test with more. 



2 METHOD 

The microlensing effect on the brightness of supemovae at redshifts 
~ 1 is not necessarily accurately modelled by a single-scattering 
event. We therefore simulate the lensing of distant supernovae with 
a ray-tracing program. We assume the Universe is populated with 
compact dark matter candidates of a single mass. The results are in- 
dependent of the mass, provided the source is small compared to the 
Einstein radius (for canonical expansion speeds, at maximum light 
this is ~ lO^^m, which requires M > 10~^Mq). The masses are 
confined to a discrete set of lens planes, and rays are passed from 
the observer to the source plane, being bent at each stage by the 
lenses in each plane. The bend angle is the vector sum of the bends 
from individual lenses, and is computed using the Bames-Hut tree 
code (Barnes & Hut 1986). The bend angle for each lens is directed 
towards each lens, with magnitude a = ^jjg^g ^ jjjg impact 
parameter. The rays are traced to the source redshift, and the solid 
angle subtended by each of the distorted pixels is computed. Since 
surface brightness is conserved by the lensing process, the ratio of 
the solid angles of the image plane pixel and distorted source plane 
pixel gives the magnification. Full details of the simulations will 
appear elsewhere (Minty et al, in preparation), but in summary we 
use 400 — 5000 lenses in total, randomly placed on 5 planes located 
to give equal numbers of lenses per plane. The random hypothesis 
is justified in Holz & Wald (1998) and Metcalf (1999). Backwards 
rays are sent out on a 1024 x 1024 grid (see Fig. hh, with rays sepa- 
rated by a small fraction of the Einstein radius of lenses on the near- 
est plane, ensuring the pattern is fully sampled. In computing the 
bend angle, the lens pattern is assumed to be periodic. We perform 
tests varying the number of sources, size of lens planes, number 
of planes, tree-code parameters, grid size, to ensure robustness of 
the results. The need for numerical simulations is illustrated in fig. 
^ which shows a typical lightcurve for lensing in an Einstein-de 
Sitter Universe, where the dark matter is all in lenses. The optical 
depth to lensing is high, and the magnification is not well described 
by a single-scattering event. We consider four separate models, one 
with no lensing, one an Einstein-de Sitter model with all matter in 
lenses, and two flat models with matter Sim = 0.3, and cosmolog- 
ical constant Qa = 0.7. These two models differ in the proportion 
of the matter in lenses. The models are detailed in table 1. We per- 
form two analyses here, a preliminary study of the 42 high-redshift 
supemovae in the Supernova Cosmology Project (Perlmutter et al. 
1999), and a second study investigating whether SNAP should be 
capable of detecting MACHOs unambiguously. In the former case, 
we take the errors from the obsevations. For the SNAP study, we 
assume that there is an intrinsic variation in supernovae properties 
of 0.157 mag (Perlmutter et al. 1999), and a measurement error of 
0.08 mag (Metcalf & Silk 1999). These errors (in magnitude) are 
assumed to be gaussian, and independent of redshift and cosmol- 
ogy, although these assumptions could be relaxed if desired. Model 



Figure 2. A sample simulated lightcurve of a supernova aX z = 1.7 in 
an Einstein-de Sitter universe dominated by MACHOs. The timescale is set 
by the mass of lenses and their velocity dispersion; for the purposes of this 
figure, these are 1O~*M0 and 300 km s~^. 



Model 


Matter density 


Cosmological constant 

Qa 


Density in lenses 


1 









2 


0.3 


0.7 


0.1 


3 


0.3 


0.7 


0.3 


4 


1.0 


0.0 


1.0 



Table 1. Models considered. Model 1 includes scatter only from measure- 
ment enor and intrinsic variations. 



1 contains only these components of variance. Fig. |3| shows the 
scatter in the magnitudes induced by lensing for models 2-A, as a 
function of redshift. Given that the intrinsic plus instrumental scat- 
ter is ~ 0.2 mag, we see that the additional variance from lensing 
becomes very significant at redshifts > 1. 

2.1 Statistics 

SNAP proposes to obtain A'^ ~ 100 — 300 supemovae per red- 
shift interval of Az = 0.1 between z = 0.4 and z = 1.2, with 
small numbers at lo wer and higher red shift (up lo z ~ 1.7; see 
the Science Case in ittp://snap.lbl.gov for full details). For statis- 



tical analysis, we use the r.m.s. in supernova magnitude in each 
redshift bin. The discriminatory power is thus determined by how 
accurately the r.m.s. can be estimated from A'^ supernovae. This 
can be calculated analytically for gaussian distributions, but the 
distribution of magnifications induced by lensing is far from gaus- 
sian (see fig. We therefore draw N magnifications at random 
from the lensing simulations, noting that the probability of a su- 
pernova lying in a distorted pixel in the source plane is propor- 
tional to the solid angle of the pixel. These magnifications are ap- 
plied to magnitudes drawn from gaussian distributions of variance 

= O-fntriusic + (^observational- ^OT the prOpOSed SNAP SUr- 

vey we take (Jintrinsic = 0.157 mag (Perlmutter et al. 1999) and 
o observational = 0.08 mag (Mctcalf & Silk 1999). For the current 
data, we take a to be the individual estimated r.m.s. (Perlmutter et 
al. 1999). We do this repeatedly, and compute the distribution of 
the sample variance as a function of redshift. Some representative 
distributions are shown in fig. ^. We compute the probability distri- 
bution of the variance estimator Di — (nij^i — rhi)'^/ {N — 1), 
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Figure 3. Mean and r.m.s. deviation of supernova brightness vs. redsliift, 
due to lensing alone. Note tliat tlie scatter has been symmetrised. The mod- 
els have, from the top, (fio, Oa, f2i,ens)=(0.3,0.7,0.1), (0.3,0.7,0.3) and 
(1.0,0.0,1.0). 




Mognificotion (mog) 



Figure 4. The distribution of magnifications arising from lensing in the 
three lensing models at redshift z = 1.7. Models 2-4, with increasing 
f2 Lens 7 from the left. 



Figure 5. The distribution of the r.m.s. scatter of N = 150 supernovae at 
redshift z = 1.7, for the four models (1^, from the left). 



where the sum extends over j = 1 . . . N supernovae in a redshift 
bin i.fhi is the mean magnitude in bin i. We do this by Monte-Carlo 
simulation of A*' supernovae drawn from ray-tracing simulations. 
This gives us, for a given model Mk, the probability of obtaining a 
set of variances {Di}, 



p{{D,}\Ah,N) =Y[p{D,\Mk,N) 



(1) 



since the bins are independent. From now on, we suppress the A'^ in 
these probabilities. 

A complication arising for our analysis of future surveys is 
the inevitable absence of data. We seek the probability of deducing 
the correct model Mcor- Marginalising over the (unknown at this 
stage) true model, this is 



k k 

(3) 



Vp(Mco,.|Mfe) 

k 



where the last step follows if we assume equal prior probabilities 
for the true model. is the number of models considered. 

To compute the conditional probability in the last equation, 
we use the distribution of sets of variances {Di}, given that the 
true model is Mt : 



p{M,„.\Mt) = / d{D,}p{M,or\{D,}))p{{D,}\Mt) 



Using Bayes' theorem, 

p{{D,}\Mk)p{Mt:) 



p(Mfe|{A}) = 



(4) 



(5) 



we see that the probability of getting the correct model, given a set 
of data is 



p{Mcor\{D,}) = 



p{{D,}\Mt)piMt) 
E,p({A}|M,.MM,) 



(6) 



where the evidence p{{Di}) cancels out top and bottom. If we as- 
sume uniform priors for the models, the probability simplifies, and 
substitution into (W) gives 



p{PIcor\Mt) = 



pi{D,}\Mt 



J2kPi{D,}\Mk)p{Mk) 



d{A} 



(V) 
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Redshift 


Number of supemovae 


Variance 


0.45-0.55 


11 


0.064 


0.55-0.65 


8 


0.170 


0.65-0.75 


3 


0.076 


0.75-0.85 


3 


0.043 



Number of supemovae per Az = 0.1 



Table 2. Estimated variances in supernova apparent magnitudes for four 
redshift bins. Data from the Supernova Cosmology Project (Perlmutter et 
al. 1999). 



Model 






^Lens 


Relative 


SN1997K 










likelihood 


removed 


1 









0.58 


1.0 


2 


0.3 


0.7 


0.1 


0.92 


0.88 


3 


0.3 


0.7 


0.3 


0.96 


0.66 


4 


1.0 


0.0 


1.0 


1.0 


0.32 



Table 3. Relative likelihood for the four models, from the high-redshift 
supemovae observed as part of the supernova cosmology project. 



In practice we approximate the integral over sets of data by a set of 
Nr random drawings (labelled by a): 



Aft 



p{{Di}^\Mt 



J2kP{{D,UM,)p{M,) 



(8) 



We compute this probability for varying numbers of supemovae per 
unit redshift interval. 



3 RESULTS 

For the existing 42 high-redshift supemovae published from the Su- 
pernova Cosmology Project (Perlmutter et al. 1999), we can com- 
pute the relative likelihood of the four models considered. For the 
data, we use the stretch luminosity-corrected effective B-band mag- 
nitude. We subtract from it the expected magnitude for model 2 (or 
3), and compute the variance (using the standard estimator) for 4 
bins between z = 0.45 and z = 0.85. The reason for the subtrac- 
tion is that the redshift bins for which we have microlensing ampli- 
fication distributions are quite broad, and the mean apparent magni- 
tude is expected to vary substantially over the bin. For the Einstein- 
de Sitter model, the additional variance from subtracting the wrong 
evolution is negligible. The variances in the four bins are shown 
in table ^. Various effects conspire to make the 0.55 < z < 0.65 
bin the crucial one: first, at lower redshift, the lensing makes little 
difference to the variance; second, at higher redshift, there are very 
few supemovae in the bins; finally, the observed variance is high in 
the second bin. Thus the second redshift bin favours models with 
more microlensing, since the expected no-lensing variance is only 
0.04. Combining the results gives the relative likelihoods shown in 
table ^, where the likelihood of the favoured model is set to unity. 
We see that the variances in the data have a slight preference for 
a population of MACHOs, but the significance is low. In fact re- 
moving a single supemova (SNI997K) from the dataset reduces 
the variance to 0.053, reversing the conclusions and making the 
no-lensing model the preferred choice (see table ^). There is little a 
priori justification for removing this point. In fact it is anomalously 
faint, so it is not a good candidate lensing event. Thus, unsurpris- 
ingly, we are unable to leam much from this test on current data. 
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0.998 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 
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True Model: 2 (Q™ = 


0.3, Ua 


— 0.7, f^Z/cns 


= 0.1) 




1 


0.004 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


2 


0.927 


0.986 


0.996 


0.998 


0.999 


1.000 


1.000 


3 


0.069 


0.014 


0.004 


0.002 


0.001 


0.000 


0.000 


4 


0.000 


0.000 


000 


0.000 


0.000 


000 


000 


True Model: 3 (Q™ = 


0.3, Ua 


— 0.7, f^Z/cns 


= 0.3) 




1 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


2 


0.065 


0.014 


0.004 


0.002 


0.001 


0.000 


0.000 


3 


0.838 


0.949 


0.977 


0.990 


0.996 


0.999 


1.000 


4 


0.097 


0.037 


0.019 


0.008 


0.002 


0.001 


0.000 


True Model: 4 (Q„ = 


1.0, f^A 


— 0, ^Lens — 


1.0) 




1 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


2 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


3 


0.088 


0.032 


0.016 


0.006 


0.003 


0.001 


0.000 


4 


0.902 


0.968 


0.984 


0.994 


0.997 


0.999 


1.000 



Table 4. Probability of deducing, from the four considered, correct and 
incorrect models given the true model. Probabilities < 5 X 10"'* are set to 
zero for this table. The redshift distribution is assumed to be uniform across 
O.K 2 < 1.7. 



Future experiments should be able to do this task well. In ta- 
ble ^ we show the conditional probabilities of selecting models, 
given a correct model, for experiments with a flat redshift distribu- 
tion of supemovae in the range 0.1 — 1.7. Probabilities less than 
0.0005 are set to zero in the table. In fig. ^ we show how the prob- 
ability of obtaining the correct model changes as we increase the 
number of sources per Az — 0.1 bin. We assume here a uniform 
prior: i.e. all models are equally likely a priori. SNAP'S expected 
redshift distribution is not uniform, but rises towards z — 1.2 with 
a small number of supemovae extending to z = 1.7. Using the ex- 
pected redshift di stribution from a o ne-year SNAP experiment (see 
Science Case in ittp://snap.lbl.gov), the conditional probabilities 
of selecting correct and incorrect models are given in table ^. Ba- 
sically SNAP should be able to do this test very well; there is only 
a very small possibility of confusing the ^Ilbus = 0.3 model with 
higher and lower lens densities. Combined with a uniform prior for 
the models, the probability of selecting the correct model is 99.9%. 
In this paper, we have only considered a discrete number of lensing 
models. A full Fisher matrix analysis of the error expected on the 
MACHO contribution is possible if many lensing simulations are 
mn, or a good fitting formula is used for the magnification proba- 
bility, but it is clear from this study that SNAP should be able to 
constrain the dark matter nature extremely well. 



4 CONCLUSIONS 

We have demonstrated how the scatter in the supernova fluxes can 
be used to support or rule out models with MACHO dark matter. 
Lensing by MACHOs increases the variance in the fluxes, and for 
models with a high density in MACHO dark matter, the variance 
can be increased by a factor greater than 3 at redshifts accessible 
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50 100 150 



Number of Supernovae per Az = 0-1 

Figure 6. The probability of obtaining the correct model, plotted against 
the number of supernovae per redshift interval h.z = 0.1. In all cases it is 
assumed that supernovae are observed up to z = 1.7. 
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True Model 
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p(Trial Model] True Model) 
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1.0 




2 


2 X IQ-" 
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0.0 
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0.0 
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0.9999 




3 


6 X 10-5 
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0.0 
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0.002 




3 


0.994 




4 


0.003 


4 


1 


0.0 




2 


0.0 




3 


0.0004 




4 


0.9996 



Table 5. Probabilities of proposed 1-year SNAP mission distinguishing be- 
tween tensing and no-lensing models. Model parameters are given in table 



by supernova searches. Simply put, if the observed variance in su- 
pernova magnitudes is too small, it can eliminate MACHOs as the 
dominant dark matter candidate. The existing data from the Super- 
nova Cosmology Project data do not rule out MACHOs, and con- 
clusions about whether MACHOs are preferred or not are sensitive 
to inclusion or exclusion of individual supernovae. Future planned 
surveys, such as the proposed SNAP survey should be able to dis- 
tinguish models readily, and in principle provide an accurate mea- 
surement of the quantity of MACHO dark matter. 
Acknowledgments 

We are grateful to Isobel Hook and Andy Taylor for useful 
discussions, and to the High Performance Computer Group, Hitachi 
Europe Ltd, on whose machines some of the simulations were run. 



© 0000 RAS, MNRAS 000, 000-000 



.J" 




z- 

(&DLJU) UU V 



o 

LTD 



O 

o 



o 



o 

<1 

ft 

a; 
cd 
> 
o 

a; 

ft 

m 

«^ 
o 

a; 



T 96'0 6'0 



